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Abstract 

Starting from the standard and commonly used relativistic nucleon wave func- 
tion, we show that residual interactions based on gluon or Goldstone boson 
exchanges generate additional components leading to configuration mixing. 
This result suggests that realistic nucleon wave functions can be expected to 
be far more complex than seemingly successful form factor fits in relativistic 
quark models would have one believe. 
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I. INTRODUCTION 

The nonrelativistic quark model (NQM) explains qualitatively many of the strong, elec- 
tromagnetic and weak interaction properties of the nucleon and other octet and decuplet 
baryons in terms of three valence quarks whose dynamics is motivated by quantum chro- 
modynamics (QCD), the gauge field theory of the strong interaction. Lacking presently 
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a complete solution of QCD at low and intermediate energies, where nonperturbative ef- 
fects are essential, quark models are useful as they provide physical insight and can easily 
be related to empirical analyses. However nonrelativistic approaches face difficulties, when 
quantitative results are required, e.g., at larger momentum transfers, or even at q 2 = for 
some weak interaction observablesQ 

A frame for a systematic approach to a relativistic few-body theory is provided by light 
front dynamics 0. This is intuitively appealing since it uses a Hamiltonian concept and 
is thus easy to interprete. As in the nonrelativistic case the phenomenon of confinement 
is modeled by a confining force that can be chosen to be three dimensional. Residual 
interactions between quarks related to "hyperfine" splittings may be calculated in light cone 
time ordered perturbation theory. By now several residual interactions have been used 
in the context of nonrelativistic models, gluon type, Goldstone boson type and instanton 
type interactions that lead to characteristic spin, isospin, and orbital mass splittings 
in the baryonic spectrum |§-f|. All these interactions induce configuration mixing that 
beside mass splitting shows qualitatively new effects. None of these interactions have been 
considered in the light front framework so far and merely "ground state" wave functions 
of the nucleon have been used. Since a complete solution of the relativistic three-quark 
problem on the light front is still missing, we provide a further step towards this aim and 
incorporate the Goldstone boson and the gluon exchange: We calculate the spin-isospin 
matrix elements that appear in the three body equation. We show that both interactions 
lead to configuration mixing. 

Due to the spontaneous chiral symmetry breakdown (%SB) of QCD, the effective de- 
grees of freedom at the scale Aqcd are quarks along with Goldstone bosons. Chiral field 
theory involves the effective strong interaction commonly used in chiral perturbation the- 



One of the early and obvious examples for a relativistic effect is the parity violating weak current 
that leads to the axial coupling of the nucleon to the W boson, g& . 
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ory (%PT 0) and applies at scales from Aqcd up to the chiral symmetry restoration scale 
~ A x = 47rf n ~ 1.17 GeV, where f n = 0.093 GeV is the pion decay constant. It provides 
a scenario inside hadrons at distances that are smaller than the confinement scale A.qq D , 
but larger than distances where perturbative QCD applies, in which the effective degrees 
of freedom of the strong interaction are dynamical quarks (with a mass depending on their 
momentum) along with the SU(3) octet of Goldstone bosons comprising the pions, kaons 
and the eta meson. Chiral field theory dissolves a dynamical quark into a current quark and 
a cloud of virtual Goldstone bosons. 



II. RELATIVISTIC NUCLEON WAVE FUNCTION 

Within the nonrelativistic quark model (NQM) baryons are described as Fock states of 
three quarks - in the simplest case. Although not derived from first principles of QCD, the 
bulk properties of the nucleon can be described in a qualitative fashion. Meanwhile experi- 
ence in quark modeling and various experiments at high energies established the necessity 
of using a relativistic description of the nucleon. The difficulties of a proper relativistic 
description of bound states are also obvious and, to date, neither a general solution exists 
nor one that would satisfy all required conditions. Among the promising candidates for a 
covariant description of a few-body bound state is one utilizing the light front form. In 
the context of the two-body problem the connection between the four dimensional covariant 
form of the Bethe-Salpeter amplitude and the three-dimensional light front wave function 
has been described explicitly and their differences have been discussed in Ref. ||. 

So far the three quark bound state problem on the light front for the nucleon has been 
tackled utilizing several (sometimes drastic) approximations, e.g. spinless quarks within a 
Weinberg-type or a Faddeev-type equation |],|1(J. To evaluate currents with appropriate 
boost properties the wave functions can be constructed by rewriting the nonrelativistic 



results in light front variables using Melosh rotated Pauli spinors [11 



A particular problem in the three quark system is the proper modeling of the confining 
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force that presently seems to require a three dimensional form. 

Lacking at present a complete solution our strategy is as follows. For the time being we 
assume that the three quark wave function can be expanded into a set of square integrable 
functions. As an example harmonic oscillator functions have widely been used in the past. 
Indeed one may assume a three dimensional harmonic oscillator confining potential U that 
will be added to the mass (squared) operator, Mq, of the free three quark system following 
the Bakamjian-Thomas construction. This leads to harmonic oscillator functions {i? n } for 
the configuration/momentum space. For the spinless case this approach has been shown to 
lead to the proper Regge trajectories, both for mesons and baryons, provided some conditions 
on the parameters are fulfilled 0. 

The proper coordinates to describe the evolution of a particle in light cone dynamics are 
the light cone time r = x° + x 3 and the coordinates (x~ , x±) = (x° — x 3 , (x%, x 2 )) arranged 
as x = (t, x , x _\_ ). The momenta of the quark in the nucleon moving frame are denoted by 
Pi and 

= (p + =p°+p 3 ,p- = m ^ + PL ,P^j ■ (1) 

In the light-cone approach all particles are on their mass shell, i.e. m 2 = p + p~ — p\ and 
the spectrum of the total momentum, P + , is bounded from below. For all particles pf > 0, 
and the kinematically invariant light cone fractions are Xi = pf /P + . Since light cone quarks 
are off the p~ shell the three-body proper coordinates are given in terms of the + and _L 
components only, 

i 

X 2 Pl ~ Xip 2 

<?3 = ; , 

Xi + x 2 

Q3 = (1 - 3:3)^3 - ar 3 (pi +p 2 ), (2) 

where the first equation is the total momentum P and the other two are coordinates for the 
pair qs and for the odd quark Q3 (given in the (12)3 coupling scheme, viz. channel 3. Other 
channels are given by cyclic permutation of 1,2,3.). Note that qf = Qf = 0. 

4 



Several boosts are kinematic, whereas rotations (L±) are not. The kinematic light cone 
generators are transitive on the null plane, so that the wave function is known everywhere, 
once it is known in the rest frame of a bound few-body system. Finally, the total momentum 
separates rigorously from the internal motion in momentum space, which is not the case in 
the instant form, where all boosts are interaction dependent and rotations are kinematic. 

By now the structure of the ground state (i.e. without additional momentum excita- 
tions and with symmetric momentum space wave function) in the light front formalism is 
well understood and tested. In the rest frame of the three particle system with free total 

o 

momentum, P— (Mo,Mo,Oj_), the nucleon ground state may be written in the following 
way: 

tfjv(123) = ^(123) R (k h± , fc+), (3) 

where 1 = {k±, X±, r 1; . . .} denotes the momentum k±, helicity Ai, isospin T\ of particle 1 and 
ipN is the relativistic wave function depending on spin and isospin. The Dirac structures 
necessary to specify the spin-isospin components are given in Table |. However, because 
of the symmetry restrictions (S3 and parity of the total wave function) the number of 
independent ground state components reduces to three linearly independent symmetric spin- 
flavor states (see e.g. [0 and refs. therein), viz. Q 

2 

<M123) = c*4(123) Ul u 2 us. (4) 

K=0 

The quark spinors are denoted by ttj = ULc(j>i, \)Xi where the light cone spinors u^ciPi, K) 
are given explicitly, e.g., in Ref. [13] and \i are the isospinors. The spin-isospin wave 



functions invariants are given by 



2 Different momentum wave functions, Rq , for various spin-flavor invariants I K are not ruled out 
and would, in fact, be needed to describe different momentum dependence of the proton's charge 
and magnetic form factors according to recent measurements at the Thomas Jefferson Research 
Laboratory (JLab) Q 
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J (123) = u llb Gul u 3 u x + (23)1 + (31)2, (5) 
h(123) = u^rGul uz lbltl ru x + (23)1 + (31)2, (6) 
J 2 (123) = ua ■ P lb Gu T 2 u 3 u x + (23)1 + (31)2, (7) 

where G = ir 2 C, and C = «7o72 is the charge conjugation matrix and T denotes transpo- 
sition. The coefficients c K have to be determined from a dynamical equation. The larger 
number of nucleonic components even in the ground state (compared to only one in the 
nonrelativistic case) is related to negative energy contributions in the wave function. This 
is similar to the deuteron as a bound two nucleon system and gives rise to P states that are 
of pure relativistic origin. In the past, the coefficients c K have been chosen in a way that 
the resulting ground state relativistic wave function is restricted to contain only positive 
energy components. As a consequence only one independent invariant spin-flavor structure 
arises. This can best been seen in the Bargmann Wigner basis [|l|], but is also clear from the 



nonrelativistic wave function that is Melosh boosted to the light cone. In its turn this partic- 
ular relativistic wave function can be said to be a straightforward generalization (based on 
Melosh rotations) of the nonrelativistic three quark s-wave function (that contains positive 
energy Fock components only). This particular linear combination, G 2 + Gq, is then given 
by 



Vo(123) = [(«i(7 ■ P + M ) lb Gul) (u 3 u N (P)) + (23)1 + (31)2 



U!U 2 U 3 , 



where M denotes the normalization and the nucleon spinor un{P) depends on total P. In 
other words, when the mixed antisymmetric nucleon spin-isospin state is Melosh boosted, 
the combination G 2 + Gg+permutations of Table | and Eq. |8] is obtained. When the mixed 
symmetric state is similarly treated another combination, G3 — G5 — Gg, is obtained. Of 
course, in the uds basis both lead to the same totally symmetric spin-isospin nucleon wave 
function. The excited states may be constructed in a similar fashion, as has been shown in 
Ref. @. 

The equation of motion for the three quark bound state ^ ^ may be written P, [r5| . |T7[ 
according to the Bakamjian-Thomas prescription (BT) as 



* N = Gf\m 2 N )(V^ + + V^)V N , (9) 

where we have assumed a two-body interaction V and utilized the three-body channel no- 
tation (i.e. = V(12) etc.). The free Green function G 3 is given by 



GT-\m 2 N ) = P + \P--Y.P 



3 

71 J 



m 



1-Z 3 2 , 1 ^2 ^- m i 



*t + ^ v^s-E— (io) 

xix 2 x 3 (l- x 3 ) p{ Xj 

that has to be evaluated at the eigenvalue m N of ^ n, and q 3 and Q 3 are four- vectors. Note 
that C?3 is independent of the specific channel chosen for the three body coordinates (here 
channel 3). The uds basis used in Eqs. [f]and || suggests a Faddeev decomposition of the wave 
function, ^jv = ^jv^ + ^iv + ^' < ' 3 ' ) shown in Fig. |l|. In this notation the Faddeev components 
of the three-body bound state are given by (a = 1, 2, 3) 

^ =GP(m 2 N )V^V N , a = 1,2,3. (11) 

A diagrammatical form of this equation is shown in Fig. [| The integrals appearing in Eqs. |9| 
and |TT] are written more explicitely (e.g. in channel 3) 

(V^ N )(1'2'3') = £ ^(l'2',12) 5 3 '3^v(123) 

123 

= / V^(x' 1 ,q' 3± ;x 1 ,q 3± ) ^(123'), (12) 

The last equation is written for P± = 0, and spin isospin indices suppressed. Note that 
x\ + x 2 + x' 3 = 1 must be fulfilled. 

This integral will now be evaluated in the following way: We introduce the nucleon wave 
function given in Eq. |3] and calculate the corresponding diagrams of Figs. using light 
front perturbation theory. 

The diagram of Fig. ||]a, e.g., represents the ir° exchange between quark 1 and 2 with 
momentum k = pi — p[ for the _L, + components. In this case the integral, Eq. [12] becomes 
(up to a Jacobian and constants), 
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f dk + d 2 k ± v^{k + ,k ± ) + 

The light front energy denominator E\i for this particular case is given by 

E 12 = m 2 N - E' 3 - E n - E[ - E 2 , (14) 

with on-shell light-cone energies Ei and E n , while the remaining vertex function denoted 
here by = v (12) will be defined in the next section. 

In the remaining part of the paper we investigate this integral and evaluate the spin 
isospin terms. Presently we focus on two types of residual interactions between the quarks, 
Goldstone boson exchange and gluon exchange. We will show that both interactions lead 
to configuration mixing in the ground state (nucleon) wave function that has not yet been 
considered for relativistic quark models. 

We note here that for the two nucleon case integrals of this type have been evaluated for 
the Bethe-Salpeter approach || and the light front approach |Tj|. The respective equations 
can be solved by iterations. Starting from a nonrelativistic wave function already the first 
order approximation generates all the relativistic components and has been quite useful to 
study qualitative features of the wave functions. The additional wave function components 
arising are indeed essential for the description of deuteron break up even at threshold energies 
as they are related to the standard relativistic pair current corrections in a nonrelativistic 
frame work, see ||. In the case of magnetic moments of the deuteron these additional wave 
function components have been shown to quantitatively agree with the standard expansion 



in terms of pair currents [19 . 



III. GOLDSTONE BOSON AND GLUON EXCHANGE 

We discuss two different types of residual interactions. One based on gluon exchange 
(vq) has been thoroughly studied in nonrelativistic models. Another based on Goldstone 
boson exchange (v 5 ) has recently been suggested as a possibility to incorporate the concept 
of chiral symmetry breaking into quark modeling |J. 
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If the chiral symmetry breakdown is based on SU(3)l x SU(3)r, then the effective 
interaction between quarks and the octet of Goldstone boson (GB) fields $j involves the 
axial vector coupling 

8 

Ant = ~ E Qd^X^q (15) 

that is well known from soft-pion physics. In Eq. [15|, the Aj, (z = 1, 2, 8) are Gell-Mann's 
SU(3) flavor matrices, and qa is the dimensionless axial vector-quark coupling constant that 
is taken to be one here. As a consequence, the polarization of quarks flips in chiral fluctu- 
ations, — > + GB, into pseudoscalar mesons of the SU(3) flavor octet of Goldstone 
bosons, but for massive quarks the non-spinflip transitions from 7±75fc ± that depend on the 
quark masses are not negligible. Let us also emphasize that, despite the nonperturbative 
nature of the chiral symmetry breakdown, the interaction between quarks and Goldstone 
bosons is small enough for a perturbative expansion to apply. 

The nucleon wave function based on the Pauli-Melosh basis contains no small Dirac 
spinor components. Hence its spin function is uniquely specified by the S-wave 3-quark 
wave function in the case of no configuration mixing. In terms of the Dirac-Melosh basis 
its spin wave function involves the combination M (j2 + Gq |T3[ of spin invariants. Many 



electroweak form factor calculations are based on this wave function. We wish to show 
first that effective interactions like pion exchange or gluon exchange generate many other 
relativistic spin invariants starting from this component alone. 

To be specific we will drop the momentum wave function of the proton and evaluate the 
term v 5 (12)ip N that is the spin-isospin part of the ground state (nucleon) wave function. 
That is, we consider a 1-2 symmetric neutral pion exchange interaction in channel 3 of the 
nucleon equation of motion (see Eq. II], Figs. and |3], tt° exch. with uud flavor part), i.e. 



u 5 (12)ViV = v s (12){if>$ + 4>ff + *PP) 

= Am 2 q wi7sWi(«275W2) («3 [P] 75 CuJ(u 2 u N ) + U 3 [P]7 5 Cu^(uiUAr)) 

= -u' 3 [P)j ■ A;Cu' 1 T (u 2 757 ■ ku N ) - u' 3 [P)j ■ Wu^{u x ^ ■ ku N ). (16) 



The positive-energy Dirac projection operator [P] = (7 • P + Mo) in Eq. ^ originates from 
the Melosh boost of the quarks to the light cone as explained in the previous section. Here, 
the isospin part is treated separately and therefore the light cone quark spinors are given by 
Ui = u(pi, Xi), v! i = u^p'i, A'j). Note that the second term in Eq. [16| is related to the first by 
1 «-> 2 symmetry. 

The instantaneous quark propagator has been included in the quark propagator (i.e. spin 



sum) by replacing the quark minus momentum component p { by the modified value |20 



p i = P — Pj — p k — k , where the bar denotes the on-shell quantity and k the pion 
momentum. For v§(12), p 2 for example, may be written as 

P2 = ~ Y.Pi) + P2 - k~ 

i 

= AP + {p 2 - p' 2 - 2k- ) + p' 2 +k~ = Ap 2 + p' 2 + fc" (17) 

AP = P--J2P~. (18) 

i=i 

Hence, when we sum over intermediate quark helicities for the first term in Eq. [Ll] using the 
kinematics p[ = p\ — k, p 2 = p 2 + k for the +, _L components in J2\ 2 u 2^2 = (7 ■P2 + r nq)/2m q 
in conjunction with the free Dirac equation for u 2 , we obtain 

2m g U2I5U2 u 2 u N = m' 2 75(7 • P2 + m q )u N 
a 2 

= m 2 757 • ku N + Ap 2 {u 2 ^ + u N ). (19) 

The small off-minus-shell second term in Eq. [19] (and elsewhere) will be neglected in the 
following for simplicity and merely indicated by ■ • •, but will be considered and numerically 
investigated elsewhere. Note that the pion momentum k is the integration variable in the 
nucleon equation of motion. 

In order to generate nucleon spin invariants in Table |, we Fierz transform each term on 
the right of Eq. [IB] using the scalar line of Table [TI|. This produces five terms which collapse 
to the following two upon using their (1 <-> 2) symmetry or antisymmetry, 
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v^\2)%l) N = ~ u x i^Cu 2 T (%[P]757 • ^7 M 7 • ku N ) 
1 



+ ^ u\a, lu Cv! 2 T (%[P] 75 7 • ka^-y ■ ku N ), (20) 
10 



whose 1-2 symmetric first matrix elements already have the canonical form of the nucleon 
spin invariants of Table |. The second part is reorganized to yield 

v 5 (12)ip N = -u'^Cu? u' 3 7 5 [ - M fcV + k 2 P^ + 2M k^ ■ k - 2k ■ Pk^UN 
+\u\a^Cv! 2 T u' 3 -f 5 \ik 2 {YP u ~ 1 U P^) + 2i{P»k v - P u k^ ■ k 



2 



+2ik ■ P{k ii 1 v - k u Y) u N (21) 



modulo small off-minus-shell terms. The first term, —Mok 2 , in the brackets gives a G3 
contribution in Eq. [21], the second and fifth give G5 and G$, respectively. Altogether we 
obtain from Eq. |2l| 

v 5 (12)ip N = -PMqu'^^Cu^ u' 3 Yl5UN ~ k 2 u[j ■ PCu' 2 T u 3 j 5 u N 

-~k 2 u 1 ia fJiU P u Cu' 2 T u'^'j^un + 2k ■ Pu'tf ■ kCu^ u' 3 ^un 
—2k ■ Pu 1 ia fMU k u Cu'2 u'^^un — 2M u 1 / y ■ kCu 2 %757 • kuN 
+2u 1 ia^P^k u Cu 2 T u 3 ^ • ku N . (22) 

In Eq. ^ the first three terms can be combined to the single proper vector-spin invariant of 
the nucleon with projector [P] from the Melosh boost for each quark, viz. 

k 2 



4M 2 



u[[P] llx C(u' 2 [P]y u' 3 [P\^ lh u N , (23) 



which has only positive energy components. The last four terms containing k dependent 
spin invariants will be addressed in the next section, where it is shown that they give rise 
to N and N*-type spin invariants through the angular k± integration. Symbolically we may 
write Eq. |22| as 

v 5 (12)(G 2 + G 6 ) ~ G 3 - G 5 - G 8 + ■ ■ ■ (24) 

in terms of the spin invariants of Table |. This says that a pseudoscalar interaction t>s(12) 
between quarks 1 and 2 acting on the nucleon spin invariants G2 and Gq is off-diagonal and 
makes transitions to the nucleon spin invariants G3, G5 and G$ so far, but more will follow. 
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We also notice from missing [P]'s next to the quark spinors in the k dependent terms 
that small Dirac components are generated by these interactions, so that the Pauli-Melosh 
basis is clearly too limited. 

The isospin matrix element corresponding to the neutral pion exchange 1*5(12) is straight- 
forward to evaluate 

(xIwaXiXxiYsXT) = !» Xi = \u), X2 = \u), xs = \d). (25) 

The same result is valid for 1 <-> 2. The isospin matrix element for the charged pion exchange 
i>5~(13) that we consider next is evaluated similarly, viz. 

-xfr-ir 2 (r + ) T Xi(xIx T ) = -2, (26) 

where = T^( T i ± "~ 2 ), t • r — r 3 r 3 — (t + t~ + t~t + ). 

To demonstrate the case for the n + exchange interaction we consider channels 1 and 
2, since in the uds basis and for the proton quarks 1 and 2 are both up quarks and no 
charged exchange is possible between them. We consider the 1-2 symmetric combination 
G {v£(23) +t>5~(32) + v£(13) +i>5~(31)] on G 2 , where the Green function symbol Qq denotes 
various energy denominators. To preserve the uud flavor order in the nucleon wave function 
we need to relabel the u and d quarks 2^3 and 1 <-> 3, respectively. Thus we start from 

\E£ + £ 3 - 2 > 5 + (23) + (E£ + E^)v+(13)] ^ N 

= Am 2 q (E~i + E~i) u' 2 ^u 2 M375M3 {u 3 [P\^ b Cu\ u 2 u N + u 3 [P]-f 5 Cu 2 WiWat) 

+Am 2 q (E~ 3 1 + E~i) ui7 5 tti u' 3 ^ b u 3 (u 3 [P]-f 5 Cu^ u 2 u N + u 3 [P}>y 5 Cu 2 UiU N ^ 
= ~( E 23 + E m) [«' 3 757 • k[P]Cu'T w' 2 7 • k^u N 

+M3757 • fc[P]7 • kCu' 2 u[un \udu) + (1 <-> 2)|di«i) 
= ~( E 23 + ^ra) [w 2 757 • k[P} lb Cu? m 3 7 • A;7 5 Mjv 

+w 2 757 • fc[P]7 • A;C-u 3 u[un \uud) + (1 <-> 2) (27) 

requiring no Fierz rearrangement for the (12)3 and (21)3 terms. The flavor dependence is 

12 



displayed on the right hand side of Eq. |27| to avoid confusion in a comparison with Eq. p~6| . 
Equation ^ will be continued in the next Section after the angular k± integration. 

IV. RELATIVISTIC COUPLED EQUATIONS 

We now turn to evaluate the integral given in Eq. [13| in first order approximation. The 
angular k± integration is performed using the identity 

J d 2 k ± a ± k ± b ± - k ± f{k 2 ) = ira ± ■ b ± J dk\k 2 J{k 2 ). (28) 

If the constant vectors a and h are the Dirac gamma matrices we write for two four-vector 
products 

/ d 2 M 7 )i • k (7)2 • kf{k 2 ) = tt( 7 )i • (7)2 / dk\k 2 f{k 2 ). (29) 

modulo 7 ± terms that we ignore at first. Here the sub- indices stand for the first or second 
matrix element in nucleon wave function components where the 7's occur. 
Using 

J d 2 k ± k ■ P 7 • kf{k 2 ) = 7 • P 7T J dk 2 ± k 2 f{k 2 ), (30) 

the remaining spin invariants in Eq. |22| are now approximated by 

J d 2 k ± k ■ P rZ' l7 • k Cul f{k 2 ) = u' l7 • PCul tt J dk 2 ± k 2 f{k 2 ), 
J d 2 k ± k ■ P u\io^k y Cul f{k 2 ) = u 1 ia llu P u Cul tt J dk 2 ± k 2 f{k 2 ), 
J d 2 k±u' x i ■ kCul U3757 ■ ku N f(k 2 ) = u'^Cu^ u z ^^u N it J dk\k 2 f(k 2 ), 
J d 2 k±u 1 ia^ u P fl k' / Cu2 W3757 • ku^ f(k 2 ) = u^a^P^Cu^ u^^un tc J dk\k 2 f(k 2 ). (31) 

At this stage we have to specify f{k 2 ) for each case in more detail. For 1*5(12) the 
kinematics p[ = p% — k, p' 2 = P2 + k mentioned above imply that E12 (see Eq. |T^p and 
Ro(q 2 , Q 2 ) will depend on k± via p' 1± ■ kj_, p 2J _ ■ kj_, and q' 31 _ ■ k± where the final momenta 
(with ') are held fixed. In fact, (Ei 2 )~ l + (E 2 i)~ 1 , Rq and M depend on q' 3 ■ k, as required 
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by translation invariance. Note that E12 contains — k± ■ p' 21 _ via E2, while E21 contains 
+k± ■ p' l± via E\, etc. which combine to q 3 ■ k± = — g 3 ■ k. 

Therefore, we may display the angular k± dependence by expanding 

R (q 2 3 , Ql) = (i2o) + q' 3 ■ k(R ) q . A + Q' 3 ■ k(R ) Qa + ■■■, (32) 

where (-Ro)o) (Ro)g 3 etc depend on the relative momentum variables q' 2 , Q' 2 and on k\ only 
without angular k± dependence. Similar expansions are valid for suitable combinations of 
energy denominators and for M (q 2 ,Ql). 

Now the last four terms in Eq. |22| in conjunction with Eq. [H] after the angular integration 
may be written as 



[ d 2 k ± k ■ Pu'tf ■ kCu'2 u' 3 j 5 u N Ro (— + — 

J V-C'12 -£^21/ 

= u' l7 • PCti? u' 3l5 u N it [ dk 2 ± k 2 (R ) (4~ + + " " " ' ( 33 ) 



f d 2 k±k ■ Pu' 1 icr^ u k v Cu 2 T u 3 j 5 ^u N Ro (— + — 

J V-C'12 -£^21 J 

= u^ia^Cu^ u' 3l ^u N n f dk 2 ± k 2 (R ) (-=?- + -=?-) +■■■, (34) 

J \£/i2 -&21/0 

/ d 2 k ± u 1 >y ■ Wu'2 M3757 • ku N R ( — + — ) 

J V-C'12 ^21/ 

= u' lltl Cu% u 3l ^u N n I dk 2 ± k 2 (R ) f J- + J-) + ■ ■ ■ (35) 



J d 2 k ± u' l ia^k v Cu' 2 T ^757 • ku N R (— + — ) 

= u\ia^Cu' 2 T u' 3l ^u N it [ dk 2 ± k 2 (R ) (J- + J-) + ■ ■ ■ . (36) 

J V-^12 -&21/0 

Before we comment on Eqs. |33] to let us continue the ir + exchange started with Eq. |27]. 
Integrating Eq. ^7] yields 

1 1 \ ( 1 1 



X2 = [ d 2 k ± L + (23) (— + — ) + ,£(13) (— + 



E-\ 3 £3 



#0^ 



-«2757m[ jP ]75C , % T u'^^un 71 f dk 2 j_k 2 (R ) (— + -j-) 

J \i^23 -£^32/0 

-u / 1 7 5 7 M [P] 75 C , u 2 T u 3 Y15Untt / dk 2 ± k 2 (R ) [ — + — 

J \£/i3 -&31/0 

-w' 2 757a[ j P]7 A Cw 3 T u'iUjvTT / dk 2 ± k 2 (Ro) f^- + 

J V-C/23 -£^32/0 
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-u[w x [Ph x Cu'f u' 2 U N 7T [ dk 2 ± k 2 (R ) fJ- + -L) 



+ 



(37) 



In Eq. |37] the first two terms can be simplified and the last two Fierz rearranged to the 
canonical (12)3 order, so that Eq. [J7|may be symbolically written as 

3„ 3 



/ d 2 k ± G (vt(23) + v£(13)) (G 2 + G e ) ~ -G x r_ + -G 3 r. 



+G 6 + ^G 7 r + + 2G 8 R+ + 



where r± , R± are radial integrals defined by 

1 



r± = it J dk 2 ± k 2 (M R ) 



1 \ / 1 1 
+ — ± — + 



E23 E32J 



-^13 -^31 / 



(38) 



(39) 



R ± = n J dk\k 2 (R \ 



23 



31/0 



(40) 



Clearly Eqs. |33] to [36] separate the last four terms of Eq. [22] into the nucleon spin invariants 
G3, G5, Gg and radial integrals. The trem with v$ now is seen to contain many other proton 
spin invariants other than G 2 - The ellipses indicate new N*-type spin invariants generated 
from higher order terms in expansions like Eq. |32], which we address next. 

When terms linear in k ■ q' 3 and k ■ Q' 3 are considered as exhibited in Eq. |32], then the 
relevant replacement of the identity Eq. ^ is 



J d 2 k ± k ■ g 37 • kf(k 2 ) = 7 ■ q' 3 vr J dk 2 ± k 2 f(k 2 ), 



(41) 



etc. This entails replacing the two 7 • k factors in expressions like Eq. [U} and [27\ resulting 
from Eq. [L3| by 7 • g 3 's or 7 • q' 3 7 • Q' 3 , etc. in the coupled radial equations of motion. The 



resulting spin invariants are characteristic of those of N*'s |]T3| with one or two quarks in 
higher orbitals. Thus, N* spin invariants in coupled equations for the nucleon originate from 
the entanglement of momentum wave functions (i. e. orbital angular momentum) with the 
3-quark spin structure. 

The first higher order terms contributing to Eq. [33], for example, are given by 



15 



9 3 ' QWil ■ p Cu'T u' 3 ry 5 u N 7T I dk 2 ± k 4 {R ) Q3 (— + — 
+Q' 3 ■ P«' X 7 ■ ggCuf u 3 ~f 5 u N n f dk 2 ^(RQ)Q z (— + — 
+q' 3 ■ Pu' xl ■ Q' 3 Cu' 2 T u' 3l5 u N 7i f dklk\R ) Q:i ( + — 



+ ■ ■ 



(42) 



'/:: 



where the second invariant vanishes when the Dirac equation is applied and the last invariant 
corresponds to a component of the negative parity N*(1535) configuration Su. However, 



upon integrating there are higher order continuations for every term in Eq. |22|. E. g. the 
second term continues with 

/ d 2 k ± k 2 R (J- + J-) = • • • + q' 3 - Q' 3 tt / dk\k\R,) Q . A (J- + J-) 

tt / d^ 4 (i?o) 93 (J- + -L) +■■•. (43) 



Clearly, these terms represent higher order contributions from the expansion of the momen- 
tum dependence of the radial wave function and energy denominators. 

Next we turn to the gluon exchange as an example of a vector interaction. Again we 
start from the 1-2 symmetric fc(12) for simplicity, 

v G (12)tp N ee Am^u'^Ui m 2 7m m 2 {u 3 [P\^Cul uxUn + u 3 [P}-y 5 CuJ u 2 u N ^j 

= -m 3 [P]7 • k'j^Cu'rf • ku N + Au 3 [P}^Cul t^Ujfp^ ■ p 2 + (1 2) 

+2u 3 [P] l5 Cu' 2 T u' l7 • p' 2 7 • ku N - 2m 3 [P] 7 • ^7 • p'i^Cu' 2 T u[u N -(1^2). (44) 

By Fierz rearrangement Eq. |34] leads to the following results with canonical (12)-spin matrix 
elements 



v G (12)ip N = -v!xJxCu2 u 3 [P\^f • ^757 7 • ku N + 2p[ ■ p 2 \2M u' 1 'y\Cv! 2 u 37 5 7 u N 
-2u' l7 • PCv! 2 u^su N + u 1 ia fjLU P u Cu' 2 u 37 5 7 M MiV 



+- u[Cu 2 u 3 [P] [ - 7 • kj • (pi + p' 2 ) 7 5 + 7s7 ■ (Pi + P 2 )7 ■ & 



+ 2 u 'ilv Cu 2 u 3 [P][-j-kj- (p[ - p' 2 )j 5 Y - 7s7 M 7 • (Pi - Pah • 



mat 



+ - u^a^Cu'T u 3 [P}[ - j ■ kj ■ (p[ - p' 2 )i^ p - 75(7^7 • (pi - p' 2 )7 • fc 



mat 
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+o u 'il5in Cu v u 3[ p ] h-ky- (pi + p' 2 )Y - ■ (pi + p' 2 )7 • k 



2 

+ 2 Us[P] [ - 7 • fc 7 • (Pi + P2) + 7 • (Pi - P' 2 )7 ■ k] un, (45) 



which contain terms linear in k and the p\. As was the case for the pion exchange, for 
gluon exchange the explicit 7 • p\ factors lead to N*-type spin invariants here also. In the 
d 2 fc_|_-integration, a la Eq. [29], the linear term 7 • k combines with q' 3 ■ k from the Rq expansion 
Eq. when 7 • (pi — p 2 ) is present as well in Eq. |45|, or with Q 3 ■ k when 7 ■ (pi + p 2 ) is 
present, etc.0 Thus, Eq. ^5] leads to the following spin matrix elements which can be readily 
manipulated into canonical form: 



u 



for the first term, i. e. the vector invariant G3 of Table [I], while the third term gives G5 and 
the fourth Gg. The second term in Eq. [15] also gives G3. For the scalar (12)-term u^Cu'^ 
the associated second spin matrix elements reduce to 

- v*[P\l ■ Q37 • (p[ + p' 2 )^u N + m 3 [P]7 5 7 • (p[ + p' 2 )>y ■ Q' 3 u N = -4M 2 (1 - x' 3 )u 3 ^ ■ Q' 3 u N , 
i.e. an N*-spin invariant. For the vector (12)-term u'^^Cu'^ they are 

-u' 3 [P]~/ ■ g 3 7 ■ (p\ - p' 2 )^Yun - u 3 [P]757 M 7 • (pi - P 2 )7 ■ % u n 

= i-±P% ■ (pi - p' 2 ) + ±q?P ■ (pi - P2)Ki 5 u N 

+AM q' 3 ■ (p[ - p' 2 )u' 3l ^u N - AM q' 3 ^u' 3l5l ■ (p[ - p' 2 )u N (47) 

i. e. (j5,(ji,(j3 N-spin and an N*-spin invariant. Terms with 7 ■ q' 3 replaced by 7 • Q 3 
are similarly reduced to the canonical forms of Table 1. The tensor (12)-term of Eq. [|6| 
generates G7 and different N*-spin invariants, the axialvector (12)-term generates Gq and 



3 Strictly speaking, a four-vector p must be replaced by p = p—Pp-P/P 2 to maintain orthogonality. 
We omit the tilde for simplicity. 
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N*-spin invariants, while the pseudoscalar (12)-term generates G 2 and three different N*- 
spin invariants. In summary the factorization into radial integrals and spin invariants for 
the gluon exchange may be written symbolically as follows 

j d 2 k ± v G (12)^ N = nj dkl(R ) (^- + ±-) Q [AM (k 2 + p[ ■ p' 2 )G 3 

-4(k 2 +p[-p 2 )G 5 + 2t p[-p' 2 G 7 



+2M q' 3 -(p' 1 -p' 2 )G 3 -2M 



- 2M 2 (1 - 4)^7 • PCu' 2 T «' 3 7 5 7 • Q' 3 u N 



2q' 3 -(p' 1 -p 2 )G 5 + 2 



(x' 2 - x'^rriq 



P-(p[-P2)Gi 



(x 2 x^r rtg _ t f . , 
1-4 



1 - x' 3 ) 
u^CuJ M3757 • (pi - p' 2 )u N 



+{Ro)q 3 + ^) [ - 2M o<^Cu' 2 T u' 3l ■ Q' 3 7 • (pi + P>iv 
+M0Q3 • (pi + p 2 )u x ^Cu 2 u' 3 u N - 2Q' 3 ■ Pu'-ffaCu? m 3 7 • (p[ + p' 2 )wjv 
+2P • (p' 2 + p 2 y i75 C< T fZ' 37 • Q' 3 u N ] +... (48) 



V. CONCLUSIONS 

The light front form provides a framework to treat few-body systems in a relativistic 
covariant way. This approach is particularly suited when boosts are important, which is the 
case, e.g., for form factor calculations. So far nucleon light front wave functions have been 
treated in a restricted way in this context, namely using positive energy projected wave 
function components only that can be directly derived from the Pauli-Melosh basis. Besides 
this restriction, the spatial part of the wave function have been assumed symmetric and in 
the "ground state" . In this paper we have addressed this issue and shown that residual in- 
teractions between the quarks lead to new components in the wave functions. These can be 
compared to configuration mixing in the nonrelativistic approach and have not been studied 
so far. They are expected to play an important role for the mass splittings, the E2/M1 ratio 
of the AN transition, the different high energy behavior of magnetic and electric form fac- 
tors of the nucleon, and other important structure information of the baryons. For example, 
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since small Dirac components are included in the spin-flavor invariants, relativistic effects 
such as some pair currents are already included in the impulse approximation. (Compare 
Refs. p|,|IB|| for the two-nucleon case.) We have investigated two currently used residual in- 



teractions, i.e. 75 coupling representing a typical Goldstone boson exchange and 7 M coupling 
representing the gluon exchange and found many new spin-isospin structures not present in 
the simple nucleon wave function. 

The uds basis used in this context suggests working in the Faddeev formalism for the 
three body system. We have evaluated a few of the integrals appearing in the respective 
equations utilizing light front time ordered perturbation theory. As a full solution of the 
problem still needs further work a next step in this direction is a one iteration approximation 
that has been proven quite useful in the two nucleon case and might be useful in the quark 
context also, because of the perturbative nature of the Goldstone boson exchange (that is 
based on chiral perturbation theory). 
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TABLES 

TABLE I. Relativistic spin invariants for the nucleon. All invariants are symmetric with 
respect to exchange of 12. With respect to spin or isospin alone G2, G4, and Gq are antisymmetric 
in the pair 12, the others are symmetric. 

Gi : 1G ® 75MA 

G 2 : 75G ® m a 

G 3 : 7^rG ® 7s7m tu a 

G 4 : 7 M 7sG <g> 7^m a 

G 5 : 7 • PtG ® 75TMA 

G 6 : 7 • P75G <g> M A 

G 7 : ct^tG ® 75<V tu a 

G 8 : ict^P u tG ® 757m^a 



TABLE II. Fierz rearrangements for scalar (S), vector (V), tensor (T), axial vector (A), and 



pseudoscalar (P) operators (matrices). For definitions see Ref. [12 |. 



S V T A P 



s 

V 
T 
A 
P 



1 


1 


8 


4 


4 


1 


1 

2 





1 
2 


-1 


3 





1 

2 





3 


1 


1 

2 





1 

2 


1 


1 


1 


1 


1 


1 
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FIGURES 




FIG. 1. Faddeev components of the wave function, see Eq. |8|. 





FIG. 3. Light cone time ordered neutral pion exchange between quarks 1 and 2. 





FIG. 5. Gluon exchange between quarks 1 and 2. 
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